Short chains at solid surfaces: wetting transition from a density functional approach 
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A microscopic density functional theory is used to investigate the adsorption of short chains on 
attractive solid surfaces. We analyze the structure of the adsorbed fluid and investigate how the 
", wetting transition changes with the change of the chain length and with relative strength of the 

("~>) ■ fluid-solid interaction. End segments adsorb preferentially in the first adsorbed layer whereas the 

concentration of the middle segments is enhanced in the second layer. We observe that the wetting 
i temperature rescaled by the bulk critical temperature decreases with an increase of the chain length. 

For longer chains this temperature reaches a plateau. For the surface critical temperature an inverse 
r-{ ' effect is observed, i.e. the surface critical temperature increases with the chain length and then 

h— -j \ attains a plateau. These findings may serve as a quick estimate of the wetting and surface critical 

1^-^ . temperatures for fluids of longer chain lengths. 

y—i \ 

r— i I. INTRODUCTION 

<4H 

Alkanes are often used as standard substances in studies of wetting and adhesion. This is connected with the indus- 
i t ' trial importance of these chemicals, being major constituents of fuels, oils and lubricants. Consequently, the literature 
aiming at theoretical description of adsorption of chain particles is very rich - see for the review Refsi^. Among 
different theoretical approaches based on application of lattice-off models, three groups of methods are of a particular 
importance. These are computer simulation techniqueSfi* 3 ^, the self-consistent field theory calculations 2, 5,6 and the 
approaches based on the applications of different versions of density functional theory (DFT)iS*2iiSiiiiiSiiiiii5iifitil. 
In recent years there have been proposed several density functional approaches for molecular systems. One of the 
q ■ earliest density functional-based studies were carried out by Woodward^ more than ten years ago. Next, Yethiraj 
and WoodwarcUSiii proposed a theory that combined a single-chain Monte Carlo simulation and a DFT scheme. 
Another class of density functional methods was based on on the interaction site formalism of Chandler et alii 2 ^. 
Finally, we should mention about theories which extends Wertheim's thermodynamic perturbation theoryi^ to the 
• case of nonuniform chain fluidsiSiii. From computational point of view the latter theories, especially the theory of Yu 
Q\ ' and WuiZ*^, seem to be particularly convenient to study the structure and thermodynamic properties of nonuniform 
. fluids composed of tangentially bonded chains. 

Majority of previous density- functional based studies of adsorption of chain particles have considered models with 
repulsive forces only and have been devoted to study of the packing of hard chains at hard walls cf. egAsLiSii&iSiiLSii. 
More recently theoretical effort has been shifted towards studies of systems with attractive interparticle and particle- 
wall forces^i* 2 ^ 3 * 2 ^^* 2 ^. Particularly worth to notice is the recent publication of Miiller et ali^iS^, who have reported 
the results of Monte Carlo, density functional and self-consistent field theory investigations of surface and interfacial 
properties of a molecular fluid composed a short linear chains (built of 10 segments). Their approach has been based on 
a model, according to which the molecules are treated as spherical sites connected by springs and with site interacting 
with other sites and with the surface via Lennard- Jones type potentials. Wertheim's thermodynamic perturbation 
theory has been applied to model connectivity of the segments. They have found 2 ^ that for the theory to be accurate, 
it is important to decompose the free energy into repulsive and attractive contributions with different approximations 
(e.g. different weighted densities) for the two parts. The developed approach has been able to reproduce the results 
of computer simulations with a good accuracy. However, it also requires much of computational effort. Consequently, 
extensive studies of surface phase transitions within the framework of the theory of Miiller et would be rather 
time-consuming. In our opinion an attractive alternative for performing such studies provides theory of Yu and 

In this work we apply the density functional theory of Yu and Wui2*i^ to study adsorption of molecules built of 
freely jointed tangent Lennard- Jones segments on a solid surface. The main purpose of our work is to investigate the 
wetting transition for short chains built of from 2 up to 10 identical segments. A mean-field approximation is applied 
to calculate the contribution to the free energy functional arising form segment-segment interaction. Therefore we do 
not expect our theoretical predictions to remain in a quantitative agreement with computer simulations. However, 
we can expect that a general picture of the wetting phenomena in such systems should be correctly predicted at a 
qualitative level. 

The evaluations of surface phase diagrams requires the knowledge of the bulk phase diagrams for the fluids in 
question. Numerous calculations of the bulk thermodynamic properties for models similar to that applied here 
have been reported in the literature, see e.gi 25 i 27 i 28 i 29 i 30 i 31 However, despite this similarity, there also exist several 
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differences, which lead to differences in the course of the liquid-vapor coexistence envelope. Therefore our surface 
calculations are preceded by the study aiming at the evaluation of the bulk phase equilibria. 

II. THEORY 

We consider a fluid of chain particles in contact with an impenetrable solid surface. Each individual particle of 
the fluid is composed of M spherical segments. Although the theory can be written down for a quite general case in 
which the chain may consist of segments of different sizes, through this work we restrict ourselves to the case of the 
segments of identical diameters a. The segments are tangentially jointed by a bonding potential acting between the 
adjacent segments of the same chain. The total bonding potential Vb(R) is a sum of bonding potentials v b between 
the segments 

M-l 

V b (R)=^v b (\r j+1 -r j \). (1) 

i=i 

where R = (ri, r2, ■ • • , ym) denotes a set of coordinates describing the segment positions. According to our treatment 
this potential satisfies the relation 

M-l 

exp[-/3VUR)] = [] 5{\v j+1 -Ti\ - a) /Ana 2 . (2) 

3=1 

Note that the model used here differs from those applied in Refs^SiSSiSS, where the so-called finitely extensible 
nonlinear elastic potential, responsible for the connectivity between the adjacent segments has been employed. In our 
approach we assume that the bonding potential, v b , is of infinitely short range. 

The fluid is in contact with a solid surface and each segment j interacts with the surface via a Lcnnard-Jones (9-3) 
potential 

v j( z j) = £ gs [( z o/zj) 9 - (za/zj) 3 ] , (3) 

where Zj is the distance between the segment j and the surface. The energy parameter, e gs , and distance parameter, 
zq, are the same for all the segments. 

Similar to previous studiesii^ we assume that the segments have a hard core forbidding any configuration that 
would lead to an overlap between two segments. An additional pair-wise attractive interaction between the segments 
is imposed. This attractive potential between two segments separated by the distance r is 

for r < er, 

»(r) = { 4e [(a/r) 12 - (<r/r) 6 ] for a < r < r cuU . (4) 
for O r cut , 

In the above r cut is the cut-off distance, and we have set its value to 3<r. The energy parameter e is independent of 
the interacting segment index. The form of the potential (@J excludes the interaction between the adjacent segments 
within the same chain. 

Within the framework of the theory of Yu and Wui£ the grand potential of the system f2 is defined as a functional 
of the local density of the fluid, p(R), 

n[p(R)} = F«[p(R)] + F ex [p(R)} + J dRp(R)(V ext (R) - m) + F att [p(R)} , (5) 

where p is the chemical potential of the fluid, F ex and are the excess free energy of hard-sphere chains and the 
ideal free energy, respectively, and F aU is the free energy due to attractive forces between particles. V ex t(R) is the 
external potential for a chain molecule assumed to be a sum of the external potential energies imposed on individual 
segments 

M 

^(R) = ^^-( ri ). (6) 
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The configurational ideal part of the free energy functional, F^ is known exactly 

md[p(R)} = h JdRp(R)V b (R) + JdRp(TL)\n(p(R)) - 1] j . (7) 
Following Yu and Wui£ we assume that F ex is a functional of average segment density defined as 

M M 

p s (r) =J2p*A*) = E dRS ( r - r M*) , (8) 



where /3 SJ (r) is the local density of the segment j of the chain. 

The free energy F ex is a volume integral over free energy density, F ex = J $(r)dr, and the free energy density can 
be split into the hard-sphere contribution, &(HS)y resulting from the hard-sphere repulsion between segments and the 
contribution due to the chain connectivity, <&(p). 

For the hard-sphere contribution we use the following expression, resulting from the fundamental measure 
theory 32 i 33 i 34 

a, l M \ , n i"2 nx ■ n 2 3 2 3 n 3 + (l-n 3 ) 2 ln(l-n 3 ) 

* ( *S) = -no ln(l - m) + — + n 2 (l - O 367rn§(1 _„ 3)2 , 0) 

where £(r) = |n 2 (r)|/n 2 (r). 

The contribution $(p) is evaluated using Wertheim's first-order perturbation theory 1 ^ 

1 - M 

*(P) = — ^— %Cln[y ( HS)(^)] , (10) 

where £ = 1 — ny 2 • ny 2 /(n 2 ) 2 , with y^ns) given by the Carnahan-Starling expression for the contact value of the 
radial distribution function of hard spheres is used 

1 - n 3 4(1 - n 3 )^ 72(1 - n 3 ) d 

with £ = 1 — ny 2 • ny 2 /(n 2 ) 2 . In Eqs. [(|§|)- (|ll[l ] the spatial dependence of all variables has been suppressed, for 
the sake of simplicity. In the above n ai a — 0,1,2,3,^1,^2 denote the weighted densities defined as the spatial 
convolutions of the average densities and the corresponding weight functions 

n a (r) = J dr'p s (r')w a (r - r') , (12) 

with the weight functions w a (r) given in Refi 3 ^. 

Finally, the mean-field attractive potential contribution to the free-energy is given by 



Fatt= E \ f drdr ' U (\ r - r '\)Psj( V )P^'( r ') 
j.j' = l.M ■' 



(13) 



If all the segments are identical, Eq. I|13|l becomes 

F att = \ J dTdv'u{\T-v'\)p s {v)p s {r'), (14) 

At equilibrium the density profile p(R) satisfies the condition 5£l{p(R)]/8p(R) = 0. This condition leads to the 
equation 

/0 (R) = exp|^-/3Vb(R)-/3 ]T X^rA , (15) 



where Xj(rj) is 



S [ F ex + F a tt ] 
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The density profile equation l|15fl can be combined with Eq. JSJ) yielding the equation for the average segment local 
density 



. M 

,(r) = exp(/3;u) / rfR^ <5(r — Vj) exp 



A I 



1=1 



(17) 



If the density distribution varies only in the z-direction the last equation can be rewritten as 

M 

p s {z) = exp(pfi)^2exp[-0X j (z)]G j (z)GM+i-j(z) , (18) 

where the propagator function Gj{z) is determined from the recurrence relation 

Gj(z) = [ dz' exph^Aj (z')} d(a ~\ Z ~ A) Gj-i (z') (19) 
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for j = 2,3,...,M and with Gi(z) = 1. 



III. RESULTS AND DISCUSSION 

We have studied adsorption of chains composed of M = 2, 4, 6, 8 and 10 segments. Prior to surface studies, we 
have evaluated the bulk phase diagrams using the theory outlined in the preceding section with the local density being 
constant and equal to the bulk density, p b . In a homogeneous system the vector weighted densities vanish, while the 
scalar weighted densities become proportional to the bulk density. From Eq. iff^l we obtain n a = £ a Mpb = £,aPs,b, 
with £3 = (cr) 3 7r/6, £2 = (c) 2 tt, £1 = cr/2 and £0 = 1. Insertion of the bulk weighted densities defined above into 
Eqs. l|9"]l-pi[l together with Eq. J7J) leads to the total free energy per unit volume, from which the pressure and the 
chemical potential can be readily obtained. The bulk phase diagrams were obtained from the condition that at the 
coexistence the chemical potentials and the pressures of a fluid in the liquid and vapor phases are equal. As pointed 
out in RefiS, Wertheim thermodynamic perturbation theory does not accurately describe the gas phase but it relatively 
well accounts for the properties of a liquid phase of provided that the system is well removed from the critical region. 

In Fig. 1 we show examples of the liquid-vapor coexistence envelopes in the reduced segment density pi — p^&cr 3 ~ 
reduced temperature T* = IcbT je plane for several values of the chain length M. Although the surface calculations 
have been carried out for M < 10, the result for M = 16 is also included here. Diamonds mark the critical points for 
all the chains lengths from M = 2 up to M = 16. For 2 << M < 10 the critical temperatures can be approximated by 
T* ,cr = 0.75764+ 0.95646 ln(M). The corresponding correlation coefficient is high and equals 0.99948. Moreover, the 
critical densities can be approximated by p c s r b = 0.26376M~ 38252 with the correlation coefficient equal -0.99955. The 
knowledge of bulk phase diagrams is necessary for the studies of adsorption from gaseous bulk phase at temperatures 
lower than the bulk critical temperature. 

We begin with the presentation of the results obtained for e gs /e — 6. For a given value of M we first solve the 
density profile equation l|18f) and then evaluate the "segment adsorption isotherm" 

POO 

T= / dz[p s (z) ~ Ps . b ] (20) 
Jo 

and the excess grand potentials AQ — £1 — f^, where is the grand potential of an uniform fluid being in equilibrium 
with the nonuniform system. The analysis of the dependence of AO on the chemical potential allows us for a precise 
location of the surface phase transitions. Note that the form of the external potential Eq.|21implies that the resulting 
wetting transition is first ordes&2ii. 

In figure 2 we show examples of the adsorption isotherms, obtained for M — 8 at two temperatures. The first 
of them, T* = 1.1 (cf. Fig. 2a) is somewhat higher than the wetting temperature. For the system in question 
the latter temperature is approximately equal to T£ w 0.845. The adsorption isotherm in Fig. 2a exhibits a well- 
pronounced prewetting jump. This jump occurs at the bulk density significantly lower than the segment density at 
bulk coexistence. Instantaneously, the adsorption before the jump is very low. Indeed, it is almost zero on the figure 
scale. 

With the temperature increase the prewetting jump becomes smaller, see Fig. 2b where we have plotted the isotherm 
evaluated at T* — 1.85. This temperature is close to the surface critical temperature, which for the investigated system 
equals to T* c w 1.965. We have also shown here the metastable branches of the isotherm. The equilibrium transition 
is marked by dashed line. 
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Figure 3 shows representative examples of the average segment density profiles, calculated for the systems from 
Fig. 2. Profiles in Fig. 3a have been evaluated at T* = 1.1 for the bulk fluid densities lower (dashed lines) and 
higher (solid line) than the density at the prewetting jump. Dotted line marks the bulk segment density for the 
corresponding liquid phase. The profile at the bulk density lower than the density at the prewetting transition is 
extremely low (the relevant curve in Fig. 3a has been multiplied by 500) and extends over approximately one layer. 
After the prewetting a thick film exhibiting a well-pronounced layered structure and extending over approximately 6 
monomer layers is developed. The film thickness is lower than the chain length and the consecutive maxima of the 
profile are at a distance approximately equal to a. Similar thick film morphology is found in simple fluids, however 
in the chain fluid considered here the first peak is relatively smaller. 

Figure 3b illustrates the film growth after the prewetting transition at T* — 1.85. Note that, contrary to Fig. 
3a, the profile before the transition (dashed line) has not been magnified. This is reflected in the smaller jump on 
the corresponding adsorption isotherm (cf. Fig. 2b). Directly after this transition, the surface film extends over 
approximately three layers. At the bulk density close to the bulk coexistence the film thickness grows rapidly, and at 
the highest density displayed in Fig. 3b spans approximately 18 monomer layers. The layered structure of the film is 
less pronounced than at T* = 1.1. Indeed, only three consecutive maxima on the profiles are observed. At distances 
from the surface larger than 4<r the thick film exhibits a plateau within which the adsorbed fluid segment density is 
close to the bulk liquid density at the bulk gas-liquid coexistence. 

Figure 4 gives a further insight into the film growth. We find that at both temperatures, T* = 1.1 (Fig. 4a) and 
T* = 1.85 (Fig. 4b) the adsorption divergence upon approaching coexistence is a power law i.e. T* oc (p^ — pl)^ 1 ^ 3 , 
where pb tC is the gas segment density at bulk coexistence. This behavior is clearly visible on the log-log plot of 
adsorption where the straight line slope equals -1/3. The same behavior has been observed for all the systems under 
study and is characteristic for the complete wetting regime in systems with long-range fluid-fluid and/or fluids-solid 
potentials^. 

Within the present approach it is possible to track down the density profiles of particular segments of a chain 
molecule. The segment local densities provide an additional insight into the structure of adsorbed fluid. In Fig. 5 
we show representative examples of the end (j = 1 or, equivalently j — 8) and middle (j = 4 or, equivalently j = 5) 
segment local densities obtained for M = 8 at T* = 1.1 (part a) and at T* = 1.85 (part b). These plots correspond 
to the systems from Fig. 2. At lower temperature (cf. Fig. 5a) the shape of the end and middle segment profiles 
is almost identical before the prewetting transition due to the fact that the density is very low (the profiles were 
multiplied by 500). After the prewetting transition we find that the consecutive maxima of both profiles are at almost 
the same distances from the surface but the first local density maximum of the end segments is higher than that of 
the middle segments. End segments seem to be more "pinned" to the surface. These findings are in accordance with 
results of Monte Carlo simulations and self consistent field theorjiS. We also note that the structure of the tails of 
the profiles is similar to that observed previously in the case of a fluid-fluid interface of a mixture of polymers built 
of tangentially jointed hard spheres^. 

Examples of the surface phase diagrams evaluated for dimers (part a) and 8-mers (part b) are plotted in Fig. 6. 
Similar to the case of simple fluids the prewetting line (solid line) joins the binodal (dashed line) tangentially in a 
wetting point (open triangle) and ends at the surface critical point (open circle). The reduced wetting temperature 
and the surface critical temperature T* c increases with an increase of the chain length. However, since the bulk 
critical temperature increases as well, it is convenient to re-scale (divide) both T* and T* c by T*. In Fig. 7 where we 
display the rescaled wetting temperature, = T^/T*, and the rescaled surface critical temperature t* c = T* c /T*, for 
M = 2, 4, 6, 8 and 10. All the calculations were carried out for e gs /e — 6. For a better visualization of the occurring 
changes we also plot here the bulk critical temperatures (marked on the right-hand side y-axis). The increase of 
the chain length leads to a decrease of t* w = T^/T* and to an increase of t* sc = T* c /T*. For M > 8 the both 
and t* c attain plateau. This information may be useful in finding a quick estimate of the wetting or surface critical 
temperatures for fluids composed of longer chain lengths. 

Finally, we address briefly the issue of the surface wettability changes with change of the strength of the segment- 
surface potential, e gs . Figure 8 presents the plot of the dependence of and T* c on e* s = e gs /e for the chains 
of the length M — 4. When e* s decreases both characteristic surface temperatures increase and the gap between 
them becomes smaller, as expected. For e* s — 2 the surface critical temperature T* c « 2.01 becomes only marginally 
higher than the wetting temperature, ps 2.03 and both these temperatures become only slightly lower than the 
corresponding bulk critical temperature. Extrapolation of the obtained results suggest that the line of the wetting 
points (triangles) would cross T* = T* for e* s « 1.8. However such decrease of e* gs would lead to a crossover from a 
complete wetting regime to a critical adsorption regime. Since in our studies we employ a simple mean-field functional 
we have not performed the relevant calculations here and leave this issue for a more elaborated approach. On the 
other hand, an increase of e gs would lead to a decrease of both T* c and temperatures. One can also expect the 
occurrence of layering-type transitions at high values of e* s and at low temperatures. We have not explored the 
layering transitions in the present work. The relevant results will be presented in a future. 
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IV. SUMMARY 



We have carried out studies of wetting behavior of fluids composed of chains built of freely jointed tangent Lennard- 
Jones monomers. The difference in the chain lengths leads to the difference in the wetting behavior. We have found 
that the wetting temperatures rescaled by the bulk critical temperature decrease with an increase of the chain length. 
For higher values of M the observed decrease becomes smaller and the curve of the chain length dependence of the 
rescaled wetting temperature, t^(M) attains a plateau for M > 9. In the case of the chain length dependence of the 
rescaled surface critical temperatures, t* c (M) an inverse effect has been observed, i.e. this temperature increases with 
the chain length and then attains a plateau for M > 9. We have also found that the thick film diverges according to 
the power law T oc (pt, c — pb)~ a ) m accordance with the general predictions for the fluids interacting with long range 
forces^. The theory allows for a detailed investigation of the adsorbed fluid structure. Our analysis has revealed that 
the end segments adsorb preferentially in the first adsorbed layer which is in qualitative agreement with the results 
of Monte Carlo simulations and self consistent field theorj^. 

It would be of interest to check our predictions against computer simulations. However such task would require 
rather long-lasting calculations, and we postpone such studies for a future. Our study cannot be considered as 
complete because we have not discussed the layering transitions transition. However, ascertaining the global surface 
phase diagrams for a system of freely jointed Lennard- Jones chains poses many difficulties and we postpone systematic 
investigations of this issue for further works. 
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Figure Captions 

Fig. 1. Bulk phase diagrams (binodals) of tangcntially jointed Lennard- Jones chains, characterized by different 
number of beads M = 2, 4, 8, 16 evaluated in the reduced bulk segment density pi - reduced temperature T* repre- 
sentation. Open diamonds indicate the locations of the critical points for chains with 2 < M < 16. The reduced units 
are denned as: p* } = p S; 6<7 3 and T* = k B T/e. 

Fig. 2. Segment adsorption isotherms of tangentally jointed Lcnnard-Jones 8-mers adsorbed on a solid planar 
surface with e* s — 6 plotted as a function of a relative undersaturation pl~pl c evaluated at two reduced temperatures: 
T* = 1.1 (part a), and T* = 1.85 (part b). Dashed line in part b marks the equilibrium prewetting transition. 

Fig. 3. (a) Average segment density profiles of tangentially jointed Lennard-Jones 8-mers before (dashed 
line, p* b = 2.0244243 x 1CT 07 ) and after (solid line, p% = 2.0244424 x 10~ 07 ) the prewetting transition evalu- 
ated for T* = 1.1. The dashed profile was multiplied by 500. (b) Average segment density profiles of tangcn- 
tially jointed Lennard-Jones 8-mers before (dashed line, pi = 0.0015) and after (solid lines, from left to right for 
p* b = 0.0016,0.0017,0.0018,0.0019,0.0020,0.0021,0.00215, and 0.002156, respectively) the prewetting transition eval- 
uated for T* = 1.85. Dotted line in both figures marks bulk segment density of the liquid phase. 

Fig. 4. Segment adsorption isotherms of tangcntially jointed Lennard-Jones 8-mers adsorbed on a solid planar 
surface with e* s = 6 (circles) plotted as a function of the logarithm of the deviation from bulk coexistence log(pj* — pi c ) 
evaluated at two reduced temperatures: T* — 1.1 (part a), and T* = 1.85 (part b). Solid line in both parts has a 
slope -1/3. 

Fig. 5. (a) Middle (solid lines) and end (dashed lines) segment density profiles of tangcntially jointed Lennard- 
Jones 8-mers before (pi = 2.0244243 x 10~ 07 ) and after (pi = 2.0244424 x 10~ 07 ) the prewetting transition evaluated 
for T* = 1.1. The profiles before the prewetting transition were multiplied by 500. (b) Middle (solid lines) and end 
(dashed lines) segment density profiles of tangentially jointed Lennard-Jones 8-mers evaluated for T* = 1.85 and for 
p* b = 0.00215. 

Fig. 6. The surface phase diagrams of tangentially jointed Lennard-Jones dimers (a) and 8-mers (b) plotted in 
the reduced bulk density - reduced temperature representation. Shown are the prewetting lines (solid lines), the gas 
branch of the corresponding bulk phase diagrams (dashed lines), the bulk critical points (filled circles), the surface 
critical points (open circles) and the wetting points (open triangles). The corresponding temperatures are marked in 
the Figure. 

Fig. 7 The dependence of the reduced bulk critical temperature T* (filled circles), the rescaled surface critical 
temperature T* c /T* (open circles), and the rescaled wetting temperature T^/T* (open triangles) on the chain length 
M. Dotted and dashed lines serve as a guide to the eye. Arrows point to the relevant temperature axis. 

Fig. 8 The dependence of the reduced surface critical temperature T* c (open circles), and the reduced wetting 
temperature T^j (open triangles) on the reduced segment-surface energy parameter e* s evaluated for the constant 
chain length M = 4. 
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